Comparing Effective-One-Body gravitational waveforms to accurate numerical data 
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We continue the program of constructing, within the Effective- One-Body (EOB) approach, high 
accuracy, faithful analytic waveforms describing the gravitational wave signal emitted by inspiralling 
and coalescing binary black holes. We present the comparable-mass version of a new, resummed 
3 PN-accurate EOB quadrupolar waveform that we recently introduced in the small-mass-ratio 
limit. We compare the phase and the amplitude of this waveform to the recently published results 
of a high-accuracy numerical simulation of 15 orbits of an inspiralling equal-mass binary black hole 
system performed by the Caltech-Cornell group. We find a remarkable agreement, both in phase and 
in amplitude, between the new EOB waveform and the published numerical data. More precisely: 
(i) in the gravitational wave (GW) frequency domain Adu < 0.08 where the phase of one of the non- 
resummed "Taylor approximant" (T4) waveform matches well with the numerical relativity one, we 
find that the EOB phase fares as well, while (ii) for higher GW frequencies, 0.08 < Mu; < 0.14, 
where the Taylor T4 approximant starts to significantly diverge from the numerical relativity phase, 
we show that the EOB phase continues to match well the numerical relativity one. We further 
propose various methods of tuning the two inspiral flexibility parameters, a^ and Wpoic, of the EOB 
waveform so as to "best fit" EOB predictions to numerical data. We find that the maximal dephasing 
between EOB and numerical relativity can then be reduced below 10~^ GW cycles over the entire 
span (30 GW cycles) of the simulation (while, without tuning them, the dephasing is < 8 x 10~^ 
cycles) . In addition, our resummed EOB amplitude agrees much better with the numerical relativity 
one than any of the previously considered non-resummed, post-Newtonian one (including a recently 
derived, non-resummed 3 PN-accurate one). We think that the present work, taken in conjunction 
with other recent works on the EOB-numerical-relativity comparison confirms the ability of the 
EOB formalism (especially in its recently improved avatars) to faithfully capture the "real" general 
relativistic waveforms. 
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PACS numbers: 04.25.Nx, 04.30.-w, 04.30.Db 



I. INTRODUCTION 



A ground-based network of interferometric gravita- 
tional wave (GW) detectors is currently taking data. Co- 
alescing black hole binaries are among the most promis- 
ing GW sources for these detectors. In order to suc- 
cessfully detect GWs from coalescing black hole binaries 
and to be able to reliably measure the source physical 
parameters, one needs to have in advance a large bank 
of ^Hemplates" that accurately represent the GW wave- 
forms emitted by these binaries. In the terminology of 
[1] one needs templates that are both effectual and faith- 
ful. The construction of faithful GW templates for co- 
alescing binaries comprising spinning black holes (with 
arbitrary masses mi, 7712 and spins Si, S2) poses a dif- 
ficult challenge. Due to the multi-dimensionality of the 
corresponding parameter space, it seems impossible for 
state-of-the-art numerical simulations to densely sample 
this parameter space. This motivates the need to develop 
analytical methods for computing (as a function of the 
physical parameters mi, m2, Si, S2) the corresponding 
waveforms. The Effective-One-Body (EOB) method [2- 
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5] was developed to analytically represent the motion of, 
and radiation from, coalescing binary black holes with 
arbitrary masses and spins. As early as 2000 [3] this 
method made several quantitative and qualitative predic- 
tions concerning the dynamics of the coalescence, and the 
corresponding waveform, notably: (i) a blurred transition 
from inspiral to a "plunge" that is just a smooth contin- 
uation of the inspiral, (ii) a sharp transition, around the 
merger of the black holes, between a continued inspiral 
and a ringdown signal, and (iii) estimates of the radiated 
energy and of the spin of the final black hole. 

The recent impressive breakthroughs in numerical rel- 
ativity (NR) [6-20] have given us access to extremely 
valuable, and reliable, information about the dynamics 
and radiation of binary black hole coalescence. It is com- 
forting (for theorists) to note that the picture which is 
emerging from the recent numerical simulations (for a 
review see [21]) broadly confirms the predictions made 
by the EOB approach. This gives us confidence in the 
soundness of the various theoretical tools and assump- 
tions used in this approach, such as the systematic use 
of resummation methods, notably Pade approximants (as 
first suggested in [1]). 

An important aspect of the EOB approach (which was 
emphasized early on [5]) is its flexibility. As was men- 
tioned in the latter reference "one can modify the ba- 
sic functions [such as A{u)] determining the EOB dy- 



namics by introducing new parameters corresponding to 
(yet) uncalculated higher PN effects. [.. .]. Therefore, 
when either higher-accuracy analytical calculations are 
performed or numerical relativity becomes able to give 
physically relevant data about the interaction of (fast- 
spinning) black holes, we expect that it will be possible to 
complete the current EOB Hamiltonian so as to incorpo- 
rate this information" . Several aspects of the EOB flexi- 
bility have been investigated early on, such as a possible 
"fitting" of a parameter (here denoted as a^), represent- 
ing unknown higher PN effects, to numerical relativity 
data [22] concerning quasi-equilibrium initial configura- 
tions [23, 24], and the extension of the EOB formalism 
by several new "flexibility parameters" [25] , and notably 
a parameter, here denoted as Wpoio, entering the Pade re- 
summation of the (energy flux and ) radiation reaction 
force. 

In view of the recent progress in numerical relativity, 
the time is ripe for tapping the information present in 
numerical data, and for using it to calibrate the vari- 
ous flexibility parameters of the EOB approach. This 
general program has been initiated in a series of recent 
papers which used 3-dimensional numerical relativity re- 
sults [26-29]. In addition, numerical simulations of test 
particles (with an added radiation reaction force) moving 
in black hole backgrounds have given an excellent (and 
well controllable) "laboratory" for learning various ways 
of improving the EOB formalism by comparing it to nu- 
merical data [30] . The latter work has introduced a new 
resummed 3 PN-accurate quadrupolar waveform which 
was shown to exhibit a remarkable agreement with "ex- 
act" waveforms (in the small mass ratio limit). In the 
present paper, we shall present the comparable-mass ver- 
sion of our new, resummed 3 PN-accurate quadrupolar 
waveform and compare it to the published results [20] 
concerning recent high-accuracy numerical simulation of 
15 orbits of an inspiralling equal-mass binary black hole 
system. We then show how the agreement between the 
two (which is quite good even without any tuning) can 
be further improved by tuning the two main EOB flexi- 
bility parameters: 05 and Wpoic- Our work will give new 
evidence for the remarkable ability of the EOB formalism 
at describing, in fine quantitative details, the waveform 
emitted by a coalescing binary. 



II. CALIBRATING i^poio, IN THE 

SMALL-MASS-RATIO CASE, FROM 

NUMERICAL DATA 



j^ <C 1, the conservative dynamics of the small mass (say 
TO2 — n) around the large one {mi ~ M) is known, 
being given by the Hamiltonian describing a test parti- 
cle /i in the background of a Schwarzschild black hole 
of mass M. On the other hand, the energy flux toward 
infinity, say F = (dE / dty^'^ , or the associated radiation 
reaction force J^rr, cannot be analytically computed in 
closed form. One must resort to black hole perturbation 
theory, whose foundations were laid down long ago by 
Regge and Wheeler [31], and by Zerilli [32] (for the non- 
spinning case considered here). The waveform emitted 
by a test particle is then computed by solving decou- 
pled partial differential equations (for each multipolarity 
(£, m) of even or odd parity w) of the form 



.2j,('r) 



d^h 



rM 



dt h 



M 






(1) 



where Vg is an effective radial potential and where the 
source term Sg^l [32-34] is linked to the dynamics-'^ of 
?7i2 — /i around toi ~ M . 

At this stage we have two options for solving Eq. (1): 
(i) use numerical methods, or (ii) use an analytical ap- 
proximation scheme for solving (1) by successive approx- 
imations. The numerical approach led, long ago, to the 
discovery of several important features of gravitational 
radiation in black hole backgrounds, such as the sharp 
transition between the plunge signal and a ringing tail 
when a particle falls into a black hole [36]. The analyt- 
ical approach to solving Eq. (1) by successive approxi- 
mations, of the post-Newtonian (PN) type, has been re- 
cently driven to unprecedented heights of sophistication 
(and iteration order). See [37] for a review. 

Our purpose in this introductory section is to illus- 
trate, on a simple case, how accurate numerical data can 
be used to optimize the resummation of PN-expanded 
analytical results. We consider the case of a particle on 
a circular orbit. The numerical solution of this prob- 
lem [38, 39] leads to an accurate knowledge of the radi- 
ated energy fiux i^ as a function of the orbital radius, or 
equivalently (and more invariantly) of the "velocity pa- 
rameter" V — {GMVlY/^. See Fig. 1 where the sofid 
("Exact") line represent the "Newton-normalized flux 
function" 



F{v) 
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') = , , ; with Fn w) = —V 
Fn[v) 5 



'\'" . (2) 



On the other hand, post-Newtonian perturbation the 
ory allows one to compute F(v) as, essentially, a Tayloi 



As a warm up towards our comparable-mass flexibil- 
ity study, let us first consider the much simpler small- 
mass-ratio case, v <^ 1. Here, v denotes the symmet- 
ric mass ratio v — mim2/{mi + 7712)^ of a binary sys- 
tem of non-spinning black holes, with masses rrii and 
1712- We also denote M = rrii + m2 ("total rest mass"), 
and fx = mim2/M ("effective mass for the relative mo- 
tion"), so that 1^ = n/M. In the small-mass-ratio limit 



^ As discussed in [30, 35] the small-mass-ratio limit of the EOB 
formaUsm leads to a generalization of the Regge- Wheeler-Zerilli 
formalism in that the dynamics of the sourcing particle fi is not 
taken to be geodesic, but is assumed to be modified by a radiation 
reaction force .?rr. The main issue of interest here is to optimize 
the resummation of the analytical approximation to .?rr , which 
is given by a badly convergent post-Newtonian expansion (known 
only to some finite order). 



series in powers of v (modulo the appearance of loga- 
rithms of V in the coefficients A„ when n > 6, except for 

n = 7) [37], say 



FTaylor(^) ^ 1 + A2V^ + AgW^ + • • • + A,,v'' + 



(3) 



It was emphasized by Poisson [39] that the successive 
Taylor approximants obtained from Eq. (3) converge 
both slowly and erratically to the numerically determined 
"exact" F{v). Subsequently, Ref. [1] pointed out that 
the resummation of the series (3) by means of succes- 
sive (near diagonal) Fade approximants led to a much 
better sequence of approximants. See Fig. 3 in Ref. [1] 
for a comparison between Taylor approximants and Padc 
approximants. The convergence of the sequence of Fade 
approximants was found to be much improved (the v^ ap- 
proximant being already very close to all its successors), 
and to be monotonic. This led to the suggestion of us- 
ing such Fade approximants also in the comparable-mass 
case, though we do not know (yet) the finite-;^ analog of 
the exact flux function F{v] v). 

The Fade resummation advocated in Ref. [1] involves 
one flexibility parameter, Wpoio, which parametrizes the 
location of the (real and positive) pole of the Fade- 
resummed F^'^'^'^{v, Wpoio) which is closest to the ori- 
gin in the complex v plane. Technically speaking, 
^Pade^^. ^p^jg) is defined as (1 — w/wpoie)~"^ times the rel- 
evant near-diagonal Fade approximant ^ of the Wpoio ~ 
modified Taylor series i^'T^y'"(w,i;poio) = FTayior(„) „ 



(^/fpole)FT-yl°'-(v) = 1 - «/Wpolc + ^2^2 + 



Ref. [1] 



advocated to use, as a fiducial value for Upoio, ''^poic = 
l/-\/3 = 0.57735 in the test-mass limit v ^ 0, and a 
slightly larger, i^-dependent value, say v^lf^{v) (moti- 
vated by Fade resumming an auxiliary "energy function" 
e(w; v) ) given in Eq. (4.8) there. Here, we point out that, 
when V ^> 0, a. slightly different choice for the numerical 
value of Vpoie can very significantly improve the closeness 
between the Fade flux and the exact (numerical) one. 

Our results are displayed in Fig 1 (a) and (b). In 
both panels, the solid line represents the "exact" result 
for the flux function F{v) as numerically computed by 
Foisson. In the upper part of Fig 1 (a) one compares 
^Exact^^^ to two different Fade (Pg ) approximants re- 
summing the same u^^-accurate (or 5.5 FN) Taylor ap- 
proximant [40]: the "standard" F^'^^iv] Wpoic = l/V^) 
and a "wpoic-flexed" [25] version of F^''^'^°{v\ Wpoio) using 
the optimized value wJj;;^*(5.5FN) = 0.5398. This choice 
of I'poio yields a Fade approximant which is amazingly 
close to the exact value. The lower panel of Fig. 1 



^ Given a certain order for tiie Taylor approximant, say i^'Taylor _ 
! + ■■■ + v^ , tile general prescription is to resum it with a near- 
diagonal Pade, P™, such that m + n = N and n = m + t with 
e = or 1. In the (exceptional) cases where such a near-diagonal 
Pade contains a "spurious pole" (i.e., a real pole between and 
""polo): one should use another choice for m and n (staying as 
close as possible to the diagonal m = n). 



(a) exhibits the differences A = F^'^'^'' - F^xact j^^. 
the two choices of Wpoio- While the standard choice of 
Wpoio (namely l/\/3 = 0.57735) leads to a rather good 
agreement (with |A| being smaller than 5 x 10^'^ up to 
V ~ 0.355, which corresponds to a radius r — 7.93GM, 
and |A| reaching 2.4 x 10^^ at the Last Stable Or- 
bit (LSO) at r = 6GM), the "flexed choice" v^^f^ = 
0.5398 ± 0.0001 yields an amazing agreement all over the 
interval < w < i;lso = l/\/6 = 0.40825 . The largest 
value of |A| over this interval is max|A| ~ 9 x 10"'', and 
is reached around v = 0.38. Note that the 4-digit accu- 
racy quoted for v^H^ = 0.5398 ± 0.0001 corresponds to 
(somewhat arbitrarily) imposing that the value of |A| at 
the LSO is smaller than about 1 x 10~^. The rounded off 
value Wpoic = 0.54 would still yield an amazing fit with 
max|A| ~ 10^3. 

In Fig. 1 (b) we explore what happens when using a 
much lower accuracy for the Taylor approximant of the 
fiux. We consider here, as an example of relevance for 
the finite i/ case, the case where one starts from a v^- 
accurate (3FN) Taylor approximant for the flux^. For 
that case the standard-choice Wpoic = l/VS still leads to a 
rather good agreement ( with |A| < 10^^ up to w ~ 0.325 
and |A|lso — 5 x 10"^), while the flexed choice Wp^j** = 
0.53 yields an excellent agreement all over the interval 
< V < fLso (with max|A| ~ 3 x 10^'^ being reached 
around v ~ 0.355). Though the closeness is less good 
than in the 5.5 FN case (3 x 10^^ versus 0.9 x 10~^), it is 
even more amazing to think that, starting from a 3FN- 
expanded flux function which (as shown, e.g., in Fig. 3 
of Ref. [1]) differs from the exact result when 0.3 ^ v ^ 
i^LSO by about 10%, a suitably flexed Fade resummation 
can decrease the difference below the 3 x 10"*^ level! 

Summarizing: In the small p limit, the value of the 
flexibility parameter Wpoio can be calibrated to yield an 
excellent agreement (from 3 x 10~^ to 0.9 x 10~^ de- 
pending on the FN accuracy) between the Fadeed flux 
function F^^'^^{v; fpoio) and the numerically determined 
"exact" flux F^^^^^{v) all over the interval < v < vlso- 
This gives an example of the use of accurate numerical 
data to calibrate a theoretical flexibility parameter enter- 
ing the EOB approach. In the following, we shall consider 
the equal-mass case, u = 1/4, and investigate to what ex- 
tent accurate numerical data [20] can be similarly used 
to calibrate the two main EOB flexibility parameters 05 
and Wpoio''- 



^ The result for the v^-accurate expansion would be very simi- 
lar and the final difference would be invisible to the naked eye. 
However, as we shall mention below, some problems with spuri- 
ous poles creep up in the near diagonal 3.5 PN Pade approximant 
when u = 1/4 and Vpoie < 0.55. Therefore we prefer to exhibit 
here the spurious-pole-free 3 PN Pade case. 

^ Note that, as already suggested in Ref. [1], one expects the "true" 
value of tipoie to depend on u. Therefore, we cannot a priori 
assume that the above best values, say ^''qj* — 0.53, will yield a 
close agreement for the flux function (or the radiation reaction) 



III. NEW, RESUMMED 3 PN-ACCURATE EOB 
INSPIRAL WAVEFORM 



After having considered the importance, for fitting 
high-accuracy numerical data, of the flexibihty parame- 
ter t;poie in the simpler small-mass-ratio limit, we wish 
to move on to the observationally urgent comparable 
mass case Av ^ 1. As we are going to see, this case 
involves two, rather than one, relevant flexibility param- 
eters: Wpoic (entering radiation reaction) and a^ (entering 
the conservative orbital dynamics). To understand the 
meaning of these parameters when Av ^ 1, let us present 
the comparable-mass version of the new, improved "ver- 
sion" of EOB which has been introduced in Ref. [30] 
and shown there to exhibit a remarkable agreement, in 
phase and in amplitude, with "exact" small mass ratio 
NR waveforms. Ref. [30] considered the small v limit, 
but with the clear methodological aim of using this limit 
to test improved EOB waveforms defined for any value 
of v. We here continue this program by comparing this 
improved EOB waveform to the recent numerical relativ- 
ity data of [20] . The improvements in the EOB approach 
introduced in Ref. [30] concern several of the separate 
"bricks" entering this approach. Indeed, it included: (i) 
a resummed, 3 PN-accurate description of the inspiral 
waveform, (ii) a better description of radiation reaction 
during the plunge, (iii) a refined analytical expression of 
the plunge waveform, and (iv) an improved treatment of 
the matching between the plunge and ring-down wave- 
forms. As the present paper will compare this improved 
EOB approach to the inspiralling NR results of [20], we 
shall only make use here of the improvement (i) . 



A. Improved, resummed 3PN-accurate waveform 

The new, resummed 3PN-accurate inspiral waveform^ 
derived in Ref. [30] takes the form (when neglecting the 
"non quasi-circular" flexibility parameters, a and 6, in- 
troduced to better represent the "plunge" which follows 
the inspiral) 



GM 22 ^ '' 



-v{r^nfF22e 
5 



-2i$ 



(4) 




(a) 




where $(i) is the EOB orbital phase, fi = $ is the EOB 
orbital frequency, r^ = rij}^'^ is a modified EOB radius ^, 
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FIG. 1: Panel (a) compares the "exact" Newton-normalized 
flux function F{v) [39] to two different Fade resummed, 
u'^^-accurate analytical flux functions: one using the stan- 
dard value Upoic = l/\/3 = 0.57735 and the other one using 
an "optimized" flexed value Wpoie = 0.5398. The bottom part 
of (a) plots the corresponding differences A = ^P^dc _^Exact 
Panel (b) plots the same quantities, except for the fact that 
it uses only w®-accurate analytical flux functions. 



in the comparable mass case v ^ f). 

Contrary to Ref. [30] where we used a ZeriUi-Moncrief nor- 
mahzed waveform 'I'22, we use here the same /i22 normahza- 
tion as Ref. [41]. They differ simply by a numerical factor: 



Rhij, 



yW- 
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-i^ 



(o)N 



The quantity r^ is such that, during adiabatic inspiral, it is re- 
lated to n by a standard Kepler-looking law H-^rJ, = 1, without 
correcting factors. However, during the plunge r^j starts signifi- 
cantly deviating from Sl~^'^ [43]. 



with tp being defined in Eq. (22) of Ref. [43], and where 
the crucial novel PN-improving factor F22 is given as the 
product of four terms 



F22it) = H,sT22f22{x{t))e 



1^22 (t) 



(5) 



Here HcS is the effective EOB Hamiltonian divided by fi 
(it describes the quasi-geodesic dynamics of the "effective 
test mass" /i), and T22 is the particularization to i = 



m = 2 of a resummed "tail correction factor" introduced 
in Ref. [30]. Its explicit expression (in the general, finite 
V case) reads 



1 



T(.-n 



r(^ + i-2ifc) ^ 



r(^ + i) 



g„feg2ifclog(2fcro) 



(6) 



where fc = GH^(,i^\mQ differs from k = mQ by a 
rescaling involving the real (rather than effective) EOB 
Hamiltonian. This "tail factor" is the exact resumma- 
tion of an infinite number of "leading logarithms" ap- 
pearing in the perturbative multipolar-post-Minkowskian 
(MPM) expansion [44-48] of "tail effects" in the {£, m) 
radiative moment. For instance, at the leading or- 
der in the monopolexmultipole interaction the radiative 
quadrupole Uij (Tr) contains a tail integral [49] 
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(7) 



while at the next to leading order it contains a tail inte- 
gral [50] 



2G''P I dTMlf\TR-T) 
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(8) 



[Here, / denotes the monopole of the source, i.e. 
/ = Madm = ^rcai]- The factor Tim resums the infinite 
series of the contributions to Ugm proportional to 



/^n jn 



drM, 



L '{rn^r) log ^— 



^ ■ (9) 



The real factor J22{x) was computed in Ref. [30] (as 
indicated in footnote 8 there) to 3 PN accuracy for all 
values of v by starting from the 3 PN-accurate multipolar 
post-Minkowskian results of Refs. [51-55]. The explicit 
form of its (PN) "Taylor" expansion reads 
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(10) 



2 log 2 + -logx. For greater 



where eulerlog(x) — 7e 

accuracy, we have added in Eq. (10) the small v limit 
of the 4PN and 5PN contributions (as deduced from the 
results of [40, 56]). 

Finally, the additional phase 522{t) is given by 

7 428 

<522 = Tjffrcal^ + -TTtT^ {H^.^f - 2AvX^'^ . (11) 

3 105 

At the time of the writing of [30] we had derived the 
full j/-dependent waveform (4)-(ll), except for the ac- 
curate value of the coefficient of vx^^"^ in Eq. (11), be- 
cause this value was not meaningful for us as it could 
be absorbed in the "non quasi-circular phase flexibility 
parameter" h included in Eq. (14) of [30]. Taking 6 = 0, 
we have recently derived from scratch the value —24 for 
this coefficient^. In the meantime, Kidder [41] has inde- 
pendently realized that it might be useful to compute the 
(, = m = 2 part of the waveform to 3 PN accuracy and 
has derived a 3 PN-accurate non-resummed (2, 2) wave- 
form [See also [42] for an earlier 2.5 PN-accurate deriva- 
tion of the (non-resummed) (2,2) waveform]. We have 
compared our result (5)-(ll) with his and found perfect 
agreement (when PN-recxpanding our result). 

Speaking of choices, there are more to be made to con- 
vert the PN-expanded amplitude factor f22^ °'' , Eq. (10), 
into a better, "resummed" EOB waveform factor. As 
said in Ref. [30], we propose to improve the conver- 
gence properties of the Taylor expansion (10) by replac- 
ing it by a suitable Pade approximant. We use the 
upper diagonal (3, 2) Pade, i.e. we use in our calcula- 



tions f22{x: v) 
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The final choice 



we need to make concern the argument x{t) of f22{x{t)) 
in Eq. (10). As discussed in Ref. [30] this argument is 
"degenerate" during the inspiral in that it can be equiva- 
lently expressed in various ways in terms of the dynamical 
variables of the system, namely (in the general finite v 
case) x — n^/^ — {r^nY = ^/r^. We emphasized that 
some choices might be better than others to automati- 
cally capture some non quasi-circular effects during the 
plunge. However, as we are here concerned with the in- 
spiral phase, we do not expect that the precise choice of 
x{t) will matter. Some preliminary checks indicate that 



^ As was emphasized in Refs. [20, 41, 57] the small additional 
phase term oc ux^''^ has anyway very little effect on observable 
quantities during the inspiral. 

* A technical remark concerning the Pade approximants we use 
(both for the flux function F(v) and the waveform /22(a;)): con- 
trary to the prescription suggested in Ref. [1], we find simpler 
(and numerically more or less equivalent) not to factor out the 
log-dependent terms (appearing at 3 PN and beyond), but in- 
stead to define the Pade approximants by considering the loga- 
rithms appearing in the PN expansion on par with the normal nu- 
merical Taylor coefficients: e.g., 1 -|- cix -I- ■ ■ ■ -I- (c^ -|- Cg log x)x^ is 
Padced by first Padeing the usual Taylor series l+cix+- ■ ■+C3X^ , 
and then replacing C3 — > c!^ -|- Cg log x in the result. 



this is indeed the case. For simphcity, we shall use here 
the argument x{t) = Q"^/^ for f22{x{t)) in Eq. (5). 



B. Effective One Body relative dynamics 

Let us briefly recall here the EOB construction of the 
relative dynamics of a two-body system (for more de- 
tails on the recent improvements in the EOB approach 
see [28, 30]). The EOB approach to the general relativis- 
tic two-body dynamics is a non-perturbatively resummed 
analytic technique which has been developed in Refs. [2- 
5, 28, 30, 43]. The EOB approach uses as basic input the 
results of PN and MPM perturbation theory, and then 
"packages" this PN-expanded information in special re- 
summed forms, which are expected to extend the validity 
of the PN results beyond their normal weak-field-slow- 
velocity regime into (part of) the strong-field-fast-motion 
regime. At the practical level, and for what concerns the 
part of the EOB approach which deals with the relative 
orbital dynamics, the method consists of two fundamen- 
tal ingredients: (i) the "real Hamiltonian" -ffroai, and (ii) 
the radiation reaction force J^^. 

It is convenient to replace the adimensionalized radial 
momentum pr (conjugate to the EOB adimensionalized 
radial coordinate r = R/M) by the conjugate Pr, to the 
"EOB tortoise radial coordinate" 



dr 



1/2 



-f 



(12) 



In terms of Pr. (and after the rescaling Hdf = Hcs/ ^i, 
Pip = P^/{^M)) the real, 3 PN- accurate Hamiltonian [4] 
reads 



HrcAr,Pr,,P<p) = M-f^roal = My 1 + 2v [HcS - 1 



with 



(13) 



Hcs{r,pr,,p^) 



\ 



pl+A{r){l + ^^ + z,$], 



El 



(14) 

where Z3 — 2z^ (4 — 3i^), and where the PN expansion of 
the crucial radial potential A{r) (= —500'^'^''™) has the 
form [2, 4] 



^Taylor (^ 



l-2u + 2iyu'^ 



^-^n^),^u* + a„.u' + Oiuu'), 



(15) 

with u = 1/r. As suggested in Ref. [5] we have 
parametrized the presence of presently uncalculated 4 PN 
(and higher) contributions to A{u) by adding a term 
+a5{h')u^, with the simple form 05(1^) ~ ac^v. [Indeed, 
it was remarkably found, both at the 1 PN, 2 PN [2], 



and the 3 PN [4] levels, that, after surprising cancella- 
tions between higher powers of u in the various contribu- 
tions to the coefficient a„(j/) = al^v + a^v^ + • • • of m" in 
A{u), only the term linear in v remained for n—2,^ 3 and 
4. Ref. [5] introduced this term with the idea that "one 
might introduce a 4 PN contribution -\-a^{v)u^ to A{u), 
as a free parameter in constructing a bank of templates, 
and wait until LIGO-VIRGO-GEO get high signal-to- 
noise ratio observations of massive coalescing binaries to 
determine its numerical value". We do not dispose yet 
of such real observations, but, as substitutes we can (as 
started in Refs. [22, 29]) try to use numerical simulations 
to determine (or at least constrain) the value of the un- 
known parameter a^. This is what we shall do below, 
where we shall compare our results to previous ones. 

As discussed in [4], the most robust choice^° for re- 
summing the Taylor-expanded function A{u) is to re- 
place it by the following Pade approximant A{u) = 
Pl[A^'^'^'^°''{u)]. Similarly, the secondary metric func- 
tion D{u) = P^[D^''^^°'' {u)] where L»Tayior(-y) ^g ^^^^^ ^^ 
Eq. (2.19) of Ref. [5]. 

The EOB equations of motion for (r, r^^,, pr^, (/3, p^p) 
are then explicitly given by Eqs. (6-11) of Ref. [28]. 
They are simply Hamilton's equations following from the 
Hamiltonian (13), except for the pip equation of motion 
which reads 



dp^ 

~dr 






(16) 



where, following Refs. [1, 3] the r.h.s. contains a re- 
summed radiation reaction force, which we shall take in 
a form recently suggested in Ref. [43] , namely 






T^ = '-^^ -^i^n^^ 



32 



./Dis(^'y; ^7 «polc) 
1 - W,p/Wpolo 



(17) 



where w^ = ^r^j, r^j = rip^'^ (with ip defined 
as in Eq. (22) of Ref. [43]). Here /dis denotes a 
suitable Pade resummation of the quantity denoted 
F'^'^y^°''(v,Vpo\c) above, i.e., the Taylor expansion of 
(1 - v/vpoic)F^''^^°'{v; v) where /'Taylor jg i-j^g Newton- 
normalized (energy or) angular momentum flux along 
circular orbits (expressed in terms of Vcuc = ^^'^ for 
comparable- mass circular orbits). Here again, to com- 
pletely deflne JF<^ we must clearly state what is the start- 
ing Taylor-expanded result that we use, and how we re- 
sum it. For greater accuracy, we are starting from 

FTayior(„. J,) = 1 + A2iu)v^ + A^iiy^ + A4{iyy 
-f A5{i^)v^ + Aq{v, logw)t;^ + A-j{v)v'^ 
-^A^{u = Q,\ogv)v^, (18) 



^ Actually, for n = 2, i.e. the 1 PN contribution to A(u), the 

cancellations even led to a complete cancellation a2(v) = 0! 
^" And the simplest one ensuring continuity with the i^ — » limit. 



where we have added to the known 3.5 PN-accurate [51- 
55] comparable-mass flux the small-mass-ratio 4 PN 
contribution [40]. Then, we use as Pade approximant 
of this (quasi-)f*-accurate expansion fuisiv, v, fpoio) = 

Pt \{l - w/fpoic)-F'^''^'°'(v; v)] ■ We indeed found that 
this specific (diagonal) Pade approximant (as well as the 
less accurate P^ one) was robust under rather large varia- 
tions of the numerical value of Wpoio (by contrast to other 
ones, such as Pj* or P^, which exhibit spurious poles 
when Wpoio becomes too small). Finally, note that, for 
integrating the EOB dynamics from some finite start- 
ing radius (or frequency) we need some appropriate ini- 
tial conditions. Refs. [3, 25, 58] indicated how to define 
some "post-adiabatic" initial conditions. In view of the 
high accuracy of the NR data of [20] (and notably of 
their extremely reduced eccentricity), we found useful to 
go beyond Refs. [3, 25, 58] and to define some, iterated 
"post-post-adiabatic" initial data allowing us to start in- 
tegration at a radius r = 15. 

Summarizing: In the comparable-mass case, the phas- 
ing, and the amplitude, of the new, resummed 3 PN- 
accurate inspiral waveform introduced in Ref [30] is given 
by inserting the solution of the EOB dynamics (given by 
Eqs. (13)-(17)) into the waveform (4)-(ll). This wave- 
form depends on two flexibility parameters, 05 and Wpoio, 
that parametrize (in an effective manner) current un- 
certainties in the EOB approach: 05 parametrizes un- 
calculated 4 PN and higher orbital effects, while Wpoic 
parametrizes uncertainties in the resummation of radia- 
tion effects (also linked to ly dependent 4 PN and higher 
radiative effects). 



IV. COMPARING THE NEW, RESUMMED 

EOB WAVEFORM TO ACCURATE NUMERICAL 

DATA 

Thanks to the recent breakthroughs in numerical rel- 
ativity one can now start to make detailed comparisons 
between EOB waveforms and numerical relativity ones. 
Working with high-accuracy data can further allow us 
to calibrate the "flexibility parameters" [25] entering ex- 
tended versions of the EOB formalism. A first step in this 
direction was recently taken by Buonanno et al. [29] by 
utilizing numerical gravitational waveforms generated in 
the merger of comparable-mass binary black holes. How- 
ever, the merger data that were used were relatively short 
(about 7 inspiralling orbits before merger when u = 1/4). 
Here, we shall instead consider the information contained 
in a recent, high accuracy, low-eccentricity simulation 
covering 15 orbits of an inspiralling equal- mass binary 
black hole [20] . [See also the previous results of the Jena 
group which covered ~ 9 inspiralling orbits [59]]. 

Boyle et al. [20] have published their results in the 
form of differences between the numerical relativity data 
and various, Taylor-type PN predictions We recall that 
there are many ways of defining some Taylor-type PN 
waveforms. Ref. [60] introduced a nomenclature which 



included three sorts of PN-based "Taylor approximants" ; 
from Taylor Tl to Taylor T3, as well as several other 
resummed approximants (Pade, and EOB). In addition, 
the definition of each such Taylor approximant makes 
two further choices: the choice of PN accuracy on the 
phasing, and the choice of PN accuracy in the amplitude. 
A fourth Taylor approximant, T4, was also considered in 
previous NR-PN comparison [26, 27, 61] and Ref. [26] had 
pointed out that it seemed to yield a phasing close to the 
NR one. Boyle et al. [20] confirmed this "experimental" 
fact and found that Taylor T4 at 3.5 PN phasing accuracy 
agreed much better with their long, accurate simulations 
than the other Taylor approximants. 

Here we shall use as EOB quadrupolar metric wave- 
form the new, resummed 3 PN-accurate inspiral wave- 
form explicitly presented above, say /if2^^(t; 05, Wpoio), 
defined in Eq. (4) and the following equations of the pre- 
vious section. Note, however, that [20] uses as basic grav- 
itational radiation variable the i — m — 2 projection of 
the corresponding Weyl curvature quantity, which is re- 
lated to the metric waveform /122 by 



/Re" 



\GM 



^22X 



it) 



dt^ 



R^ 
GM ^^ 



(t) = Ax{t)e 



^i0x(t) 



(19) 

Here, A and (/) denote the amplitude and phase of the 
curvature wave considered by Boyle et al. [20]. We have 
introduced a label X which will take, for us, three values: 
X=NR denotes the numerical relativity result of Ref. [20] , 
X=EOB denotes the improved EOB waveform presented 
in the previous section and X=:T4 will denote the "Taylor 
T4" waveform highlighted in Ref. [20] (as well as in previ- 
ous PN-NR comparisons [27, 61]) as giving a particularly 
good fit. 

The precise definition of this T4 waveform (as used in 
Ref. [20]) is as follows. As indicated in Eq. (19), ^f'^'^ is 
obtained by taking two time derivatives of a correspond- 
ing metric T4 waveform hJi- The latter waveform is 
defined (if we understand correctly the combined state- 
ments of Refs. [20] and [41]) by the following procedure. 
First, one defines a certain "T4 orbital phase" <i>T4(i) by 
integrating the ODEs 



d$ 



T4 



,3/2 



dt CM 
dx 64t^ 

'di ~ 5GM 



5 Taylor 
^ '^3.5 



(x) 



(20) 
(21) 



where 



a 



Taylor f^-^ _ 
3.5 



(x) = 1 + a2{i')x + a3{h')x ' + a4{i')x 

+ a5{iy)x^/^ + aQ{v, logx)x^ + a^{v)x''/'^ , 

(22) 



is the 3.5 PN Taylor approximant (for finite v) to the 
Newton- normalized ratio (flux-function)/ (derivative of 
energy i\xact\on)~F{v)/E'{v), which enters the adiabatic 
evolution of the orbital phase (see, e.g. [1, 60]). In 
the relevant case v = 1/4, 03 5^ ""^(x) is explicitly given 



as the quantity within curly braces on the r.h.s. of 
Eq. (45) of [20]. Having in hand the result, ^Tiit), 
nT4(i) = d<^T4:/dt = x^^/GM, of integrating Eqs. (20)- 
(21) one then defines a nPN-accurate, T4 (^ = 2,?7i = 2) 
waveform by truncating to order a;!^^ included the Taylor 
series 



gm' 



,T4 



(t) 



^g-2i*T4(t)^^^ 



- ~" 3 /2 

1 + /l2XT4 + h^Xr^^ 



n^X'^A 



u 5/2 , r 3 



T4 



(23) 



where the coefficients hn are obtained from the coeffi- 
cients hn in the Taylor-expanded 3 PN-accurate (2,2) 
waveform derived by Kidder [41] by setting to zero all 
the terms proportional to log(a;/a:;o) or log (x/xo) (but 
keeping the separate logx term entering /ig), and replac- 
ing V = 1/4. 

Note that our resummcd waveform (4)-(ll) differs from 
the 3 PN-accurate version of (23) in several ways: (i) the 
"orbital phase" evolution is given, for us, by the EOB 
resummed dynamics, (ii) we do not neglect the phase 
terms linked to log(a:;/a;o), but include them either in the 
resummed tail factor T22 or in ^22, and (iii) we resum the 
amplitude of /122 by factoring both iJeff and the modulus 
of r22, and by Padeing /22. 

Ref. [20] presents in their Fig. 19 the differences, in 
phase and amplitude (of the radiative Weyl-curvature 
components '^4), between Taylor T4 3.5/2.5 (i.e., 3.5 PN 
in phase and 2.5 PN in amplitude) and the (unpublished) 
corresponding Caltech-Cornell NR data, say 



(A0)t4NR = 014 (i) - '/'Nr(0 : 

AA\ ^ AT4(i) - ANR(i) 



A / 



T4NR 



Am{t) 



(24) 

(25) 



as plotted in the two panels of Fig. 19 of Ref. [20]. Our 
main aim here is to compare the "numerical data" (24)- 
(25) to the corresponding theoretical predictions made 
by the EOB formalism, say 



(A(/))t4EOB = '/'T4(i) - '/'EOb(0 : 

A^\ _ A^ijt) - A^o^jt) 



T4EOB 



EOB 



it) 



(26) 
(27) 



As just said, as the most complete, and best plotted data, 
concern Taylor T4 3.5/2.5, we shall use the 3.5 PN ac- 
curate Eq. (22) but only consider the 2.5 PN truncation 
of the Taylor-expanded waveform (23) (i.e. a waveform 
essentially contained in the 2.5 PN h+ and ft-x results of 
Arun et al. [62]). [As we use the T4 waveform only as 
an intermediary between the NR and EOB results, we 
are allowed to use any convenient "go between" , even 
if its PN accuracy differs from the (formal) one of our 
resummed EOB waveform]. 

To effect the comparison between NR and EOB, i.e., 
to compute the crucial difference (/)eob — 0nr, we needed 



to extract actual numerical data from Fig. 19 of [20]. We 
did that in several ways. First, we measured (with milli- 
nietric accuracy; on an A3-size version of the left panel of 
Fig. 19) sufficiently many points on the solid upper curve 
(Taylor T4 3.5/2.5 matched at Mloa = 0.1 ) ^^ to be 
able to replot (after "splining" our measured points) this 
upper curve with good (visual) accuracy. We could then 
use this splined version of eleven (adequately distributed) 
selected points on the upper curve on the left panel of 
Fig. 19 as our basic (approximate) "numerical data" . It 
gives us a (continuous) approximation to the a;4-matched 
phase difference A"^ (/)t4nr (?"''' ) = '^t4(*T4) J '/'JVfl (*"")■ 
[Here, i'^* denotes the time shifted so that i"* = cor- 
responds to the oji matching point between T4 and NR] . 
As we have separately computed 0x4(^X4) by integrating 
Eqs. (20)-(21) (and that it is easy to shift it to obtain 
'?^T4(^T4)) ^^ have thereby obtained an approximation to 
0nr(^'^'')- We then shift again the time argument to our 
basic EOB dynamical time so as to obtain (/^nrX^eob) 
(with the condition that tiJEOB(^EOB) ^ ^4 corresponds 
top" =0). 

There are then several ways of comparing (/inrX^eob) 
to the EOB phasing (/'eob(^eob), obtained from the pro- 
cedure explicated above. We wish to emphasize here that 
there is a useful way of dealing with the information con- 
tained in (any) phasing function (t>x{tx) where X=NR, 
EOB, T4, etc. Indeed, a technical problem concerning 
any such phasing function is the presence of two shift am- 
biguities: a possible arbitrary shift in (j)x {4>x — *■ (^x + cx) 
and a possible arbitrary shift in the time variable ix 
(tx ^ ix + Tx)- Similarly to what one does in Eu- 
clidean plane geometry where one can replace the Carte- 
sian equation of a curve y = y{x) by its intrinsic equation 
K = K{s) (where K is the curvature and s the proper 
length), we can here (in presence of a different symme- 
try group) replace the shift-dependent phasing function 
4>x(tx) by the shift-independent intrinsic phase evolution 
equation: dwx/dtx — ax('j^x), where ujx = d(j)x/dtx (for 
simplicity we shall use here M = mi +m2 = 1). It is also 
convenient to factor out of the phase acceleration ax{^) 
its "Newton" approximation^^, 

aN{u;) = cW'/'; c, = ^2^/3^ (28) 

5 

and to consider the reduced phase acceleration function 



^^ Ref. [20] computes various differences A^"^<p(t) = </'^2'(*T4) ~ 
<^NR,(t) where, given a "matching" frequency uim, 4'T'i i^Tl,) '^^~ 
notes a version of (^T4(tT4) which is shifted in <p and in t so 
that = A'^'n(/){t) = dA""^<j>{t)/dt at the moment f^"^ where 
d'P'NR/dt = LOm- We shall denote the four matching frequen- 
cies used in [20] as MuJi = 0.04, MLd2 = 0.05, MuJs = 0.063, 
McJ4 = 0.1 

^^ Note that we arc dealing here with the gravitational (curva- 
ture) wave frequency lu. If we were dealing with the orbital 
frequency Q = d^/dt, we would instead consider the following 
reduced phase acceleration {dQ/dt)/{Ci^i^^^''^) = Aii(fl) with 
C = {96/5)u. 
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FIG. 2: Reduced phase-acceleration curves as defined in 
Eq. (29) (with M = mi + m2 = 1). The inset highlights 
how the "tuned" EOB curve nearly coincides (for u < 0.08) 
with the NR (and T4) curves, while the "non-tuned" EOB 
one lies slightly below. 



x^ defined by 



wx 



11/3 



r(^x) 



(29) 



Independently of the label X the function a^ — > 1 as 
a; — !■ 0. The phasing comparisons then boil down to com- 
paring the various (reduced) phase acceleration functions 



-,NR 



i^), 



,T4/ 



(lu) and a^^ {uj; a5,Wpoie)- The last two 



functions can be straightforwardly (numerically) com- 
puted from the formulas written above. As for a^^'{Lu) it 
is, in principle, also straightfiorwardly computable from 
'/'nrI^nr)- However, as we do not have access to an ac- 
curate estimate of </'nr(^nr) but only to a rather rough 
cubic spline approximation to it, we can only compute 
an even rougher estimate of a^^iuj). [As is well known, 
taking (two!) derivatives of approximate results consid- 
erably degrades the accuracy.] Still, as we think that this 
is the conceptually clearest way of presenting the com- 
parison, we used the data we had in hand to compute the 
various "acceleration curves" presented in Fig. 2. 

Actually, it is instructive to include further phase- 
accelerations in the comparison. In Fig. 2 we show 
the following phase-acceleration functions (versus w, i.e., 
Mw): (i) Taylor T4 3.5/2.5, (ii) NR, (iii) a standard 
"non-tuned" EOB with 05 = (i.e. essentially the 
3 FN approximation) and Upoic = ■y?c)ic(^) [1] ^'^'' ^ = 
1/4, which corresponds to our current knowledge, (iv) 
a "tuned" EOB with 05 = 40 and Wpoio = 0.5074 (see 
below), and finally (v) the adiabatic EOB for as = 40 
and t^poie = 0.5074. Here, the adiabatic approximation 
to tti^ is that defined by the usual adiabatic approach 



to inspiral phasing (see e.g. [1]), leading to a„(cj) = 
F{v)/E'{v) where v = (w/2)^/'^ and where F is the 
Newton-normalized circular flux and E' the Newton- 
normalized derivative of the circular energy function. 
[When applying these general concepts to the EOB we 
need, as discussed in Ref. [3], to use the analytical, adi- 
abatic approximation to EOB inspiral, with, notably, 
p^ = j adiabatic ^.y^ obtained by solving dH^oB/dr = 
with pr ^ 0]. 

Several preliminary conclusions can be read off Fig. 2. 

(i) In the frequency domain (say Mlu < 0.08) where the 
T4 3.5 phasing matches well with the NR phasing, both 



the standard "non-tuned" EOB (05 = 0, Wpoie 



_ ,,D1S 
^polo 



and "tuned" EOB phasing (defined above) match well 
with the NR phasing. However, a closer look at the ac- 
celeration curves (see inset) shows that the "tuned" EOB 
phasing agrees better with NR (and T4). [The "non- 
tuned" a^°^ is slightly below a^^ (and a^"^) by roughly 
1.5 X 10-3 when Muj ~ 0.06 ]. 

(ii) For higher frequencies (0.08 < Mlu < 0.14), Tay- 
lor T4 3.5/2.5 starts to significantly diverge from the 
NR phasing. ^^ By contrast, both the standard "non- 
tuned" EOB phasing and the "tuned" EOB one continue 
to match quite well the NR phasing. This will be shown 
below by using other diagnostics than the acceleration 
curves. Indeed, when Afcj > 0.08 our "NR accelera- 
tion curve" exhibits fake oscillations which come from 
our use of a coarse approximation to NR data. The visi- 
ble "kinks" in our NR acceleration curve are due to our 
taking (numerical) second derivative of a cubic spline in- 
terpolant of approximate NR data points. We expect that 
the exact "NR acceleration curve" (computed with accu- 
rate numerical data instead of our approximate ones) will 
be a smooth curve lying close to the two EOB curves in 
Fig. 2. 

(iii) The fact that the adiabatic EOB curve diverges 
quite early, and upwards, from the full EOB curve is a 
confirmation of the conclusion derived in Ref. [3] (see 
Figs. 4 and 5 there), namely that , " in the equal mass 
case h' — 1/4 the adiabatic approximation starts to sig- 
nificantly deviate from the exact evolution quite before 
one reaches the LSO" . This further confirms the sugges- 
tion of [20] that the good early {Mcu < 0.08) agreement 
between T4 and NR is coincidental. 

Because of our lack of an accurate knowledge of a^^, 
we cannot use the acceleration curves of Fig. 2 to make 
any accurate comparison between EOB and NR data. 
In the following we shall use other tools for doing this 
comparison and, in particular, for constraining the values 
of 05 and Upoie- 



^^ Though Ref. [20] tends to mainly emphasize how well Tay- 
lor T4 3.5/2.5 agrees with the NR phasing one should note that 
the high curvature of the upper 0^4 curves when Mcu > 0.08 in 
the left panel of Fig. 19, and the subsequent fast rise of all the 
A""*(j!), are clear signals that Taylor T4 3.5/2.5 starts to signifi- 
cantly {and increasingly) diverge from the NR phasing. 
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FIG. 3: Correlation between Wpoic and as (top panel) ob- 
tained by imposing the constraint (35). The numerical ac- 
curacy with which Eq. (35) is satisfied is displayed in the 
left-bottom panel. The right-bottom panel displays the ex- 
tent to which, as a^ varies, the other ratios p^^, Eq. (34), 
approximate unity. 



The first tool we shall use consists in selecting among 
our eleven approximate points on the A'^'*0t4NR curve 
two special ones, namely 



A"^(^T4NR(i"^ 



T4NRi,''NR 

CJ4 



1809M) 



rbwd 



A"*0T4NR(iNR + 44.12M) = S. 



fwd 
4 



i 0.055 
0.01 , 



(30) 

(31) 



to which we shall refer as the (main) "backward" and 
"forward" W4 data. In addition, we also measured a cou- 
ple of selected points on the L02- and ws-matched lower 
A(/) curves. Namely, 



A" 



^-^(^NR 



A"^0(tN^R 



lOOOAf ) EE ^2 ^ -0.01 , (32) 

- lOOOM EE (^3 ~ -4.3 X 10^3 . (33) 



We can then use, in a numerically convenient way, these 
data to quantitatively compare (with an hopefully rea- 
sonable numerical accuracy) NR to EOB by considering 
four ratios, p^^, p,^^, p]^'^, pi^^ (where we recall that 



0.05, W3 = 0.063 and uja = 0.1), with 



/3w„(a5,1'polc) = 



A'^ 



yT4EOB ^t-NR 



V NR "*" ^'^rn) 



(34) 



and I 



~ UJ2, W3 and IJJ4 or uj. 



,fwd 



If our approximate measures (given in Eqs. (30)- 
(33)) of the various (5m 's were accurate, a per- 
fect match between NR and EOB would corre- 
spond to having all those ratios equal to unity: 

Pi^Mb.Vvolc) = 1, PuAo-^ ^'^pole) = 1, p^^^{a^,Vpo\e) = 1, 

and P^J_^ {0-5 T v-poia) — 1. This would give four equations 



FIG. 4: The Loo norm of the phase difference between EOB 
(when Upoie is correlated to as as in Fig. 3) and numerical 
relativity, as defined by Eq. (36). 



for two unknowns (05 and t'poic)- Even if we had ex- 
act values for the various (5„i's, we do not, however, ex- 
pect that there would exist special values of 05 and Wpoio 
for which all these ratios would be equal to one. In- 
deed, 05 and Upoic are only "effective" parameters that 
are intended to approximately mimic an infinite number 
of higher i/-dependent, resummed PN-effects. The best 
we can hope for is to find values of 05 and Vpoic allow- 
ing one to give a good overall match between (/)nr(0 and 
</'eob(0 (or "w^ and a^°^). To investigate this issue, 
it is then convenient to focus first on only one compari- 
son observable. We choose p^'^ because it is, among the 
data which we could measure with reasonable accuracy, 
the one which has the largest "lever arm". [Indeed, it 
corresponds to some weighted integral of the difference 



^EOB 



,NR 



over a significantly extended frequency in- 



terval]. Imposing the constraint 



P^f(a5,«poic) = l 



(35) 



then gives a precise way of exploring which extended 
EOB models best match the NR phasing . Note first that 
this equation could have no solutions. [For instance, if 
we were using the adiabatic approximation to EOB there 
would be no solutions]. To admit solutions is already a 
sign that EOB can provide a much better match to NR 
than T4. Then, the solutions could exist only if both 05 
and Wpoio are close to some "preferred" values. Actually, 
we found that Eq. (35) defines a continuous curve in the 
(a5,«poio) plane^**. 



Consistently with what was found for lower approximations, and 
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For all values of 05 > 0, we (numerically) found 
a unique value of Wpoie satisfying the constraint (35). 
This continuous curve is plotted in the upper panel of 
Fig. 3. When remembering that Eq. (30) is only ap- 
proximate,^^ we have to mentally replace the continuous 
curve in the upper panel of Fig. 3 by a narrow valley of 
"best fitting" values of (a5,Wpoic)- Let us first remark 
that this valley extends only on a rather small range 
of values of Upoic, around 0.55. It is comforting that 
this range includes the values that were previously sug- 
gested for VpoW- namely v'pfii^ = 0) = l/\/3 = 0.57735, 



DIS , 
'^polo 



{1^ = 1/4) ~ 0.6907, v°°Jl{iy = 0) ~ 0.54 (discussed 
above). 

To go beyond this result and see whether the other 
measurements constrain the value of 05 , we plot on the 
lower, right panel of Fig. 3 the values of the ratios Puj2, Pu^ 
and p^J^ along the p^™^ = 1 curve. As, along this curve, 
iipoic is a function of 05, the above three ratios depend 
only on 05. Ideally, we would like to find values of a^ for 
which the remaining ratios are all close to unity. [Given 
the coarse nature of our measurements, we cannot expect 
to get exactly unity]. We see on Fig. 3 that the ratio 
p^'^ is reasonably close to unity for most values of a^. 
By contrast, the two other ratios p^^ and p^j^ happen to 
have the wrong sign. This negative sign means, in terms 
of the phase- acceleration curves of Fig. 2, that around 
frequencies 102 and 0^3, aj^^(ti;) is slightly above a^'^{Lo)^ 
while it seems that aj^'-'^(tj) tends to be generally slightly 
helow a^'^{Lo). On the other hand, for larger frequencies, 
it seems clear that a^^ crosses a^* to become helow a^"*, 
and to become in rather good agreement with a^ . At 
this stage, the best we can do is to say that an overall 
best match between FOB and NR will be obtained when 
05 belongs to a rather large interval (say 10 ^ 05 < 80) 
centered around 05 ~ 40, where p^^^ and p^g are negative, 
but rather small (say —0.5 <■ p^^ <■ 0) 

To get another, potentially better measure of the 
"closeness" between NR and FOB we looked at the "£00" 
distance between the two functions 0nr(O s-'^d (/>eob(0 
on the time interval (in FOB time) 900M < iEOB < 
3460A/ (which roughly corresponds to the time interval 
plotted in Fig. 19 of Ref. [20]). More precisely, we com- 
puted the quantity 



-^NR (^ 



NRJ 



^00(05) = sup |0EOB(iEOB) 

900Af<tEOB<3460Af 

(36) 
where 4i^\^ (^nr) ^^ matched to the FOB phase at loa, 
and where the FOB was constrained to lie along the 
curve Wpoic(a5) plotted in Fig. 3 (i.e., satisfying Eq. (35)). 



for the presently computable contributions to as [25] , we expect 
that as > 0, and we shall therefore only work in the correspond- 
ing half plane. 
^^ We estimate the accuracy of our measurement result Eq. (30) to 
be such that the "backward time-shift", corresponding to a r.h.s. 
exactly equal to 0.055, is (-1809 ± 15)M. 




0.04- 

3"?°-°^ -^^.0. v' =0.5074 

-©- »0)3, ^ 

0[--A *T4E0B ♦- 

A°>*(t) 

^T4E0B 



M. 







FIG. 5; The upper panel compares various phase differences 
A"^"(jf)T4x versus time (with M = 1), Um denoting a matching 
frequency and the label X being either EOB or NR. The lower 
panel exhibits the a;4-matched phase difference between EOB 
and NR. The flexibility parameters of EOB have been tuned 
here to as = 40 and Upoic = 0.5074. 



We show this Loo norm in Fig. 4. This Figure displays 
the remarkable agreement between EOB and NR phas- 
ing over an interval where T4 exhibits a clear dephasing 
with respect to NR. Indeed, Fig. 19 of [20] shows that 
on this interval all Taylor T4 3.5 templates dephase by 
~ 0.08 radians (because of the divergence at the end, cor- 
responding to the divergence of the acceleration curves 
in Fig. 2 when lo ^ 0.08). By contrast, the dephasing 
between EOB and NR can be as small as 0.006 radians if 
30 < 05 < 52, or 0.008 radians if 10 < 05 < 80. Again, we 
find that a largish interval of 05 values centered around 
05 ^ 40 seems to be preferred (when Wpoio is correlated 
to as via the curve of Fig. 3) to give the best possible 
overall match between EOB and NR. 

To give a better feeling of how well EOB matches NR 
phasing all over the time interval explored by the sim- 
ulation of Ref. [20], we plot in Fig. 5 the superposition 
of the upper curve in Fig. 19 of [20] (i.e., the difference 
A'^'*0T4NR, as measured and splined by us) with the cor- 
responding EOB difference A'^^ (/>t4eob , for the values 
05 = 40, Wpoio — 0.5074 approximately corresponding 
to the smallest Lq^ norm in Fig. 4. We also plot the 
LJ2- and wa-matched phase differences A"^(/)t4eob and 
A"^(/)T4EOB- Apart from the slightly wrong curvatures 
of the i02- and W3- curves (for oj < 0.08), this Figure ex- 
hibits a truly remarkable visual agreement with the left 
panel of Fig. 19 of [20]. It exhibits again two facts: (i) 
the EOB phasing agrees extremely well with the NR one 
on the full time interval (900M < isOB < 3460M), (ii) 
by contrast Taylor T4 3.5/2.5 starts diverging from EOB 
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FIG. 6: Comparison between the standard, "non-tuned" EOB 
(as = 0, -Upoio = UpofoCi^ = 1/4) = 0.6907) and NR. The top 
panel shows that the gravitational wave phases 0eob and cj!>nr 
(versus time) are nearly indistinguishable to the naked eye. 
The bottom panel quantifies the small difference between the 
two. 



FIG. 7: Comparison between relative amplitude differences 
(A"^'"A/j4)t4x versus time, LUm denoting the matching fre- 
quency and the label X being either EOB (for as = 40, 
up„ie = 0.5074) or NR. 



when uj > 0.08 in precisely the same way that it diverges 
from NR. In the bottom panel of Fig. 5 we give a pre- 
cise quantitative measure of the difference between EOB 
and NR phasings by plotting the cj4-matched difference 
^EOB (i) — (Anr (^4)- This phase difference vanishes both 
when weob (^eob) ~ ^^ O^Y construction), and at the 
time t^Q = ig^B - 1809Af (by our optimized choice of 
the link fpoic = t^poio(a5), such that Eq. (35) holds). We 
see how, indeed (in agreement with Fig. 4) the dephas- 
ing remains smaller, in absolute value, than about 0.006 
radians, i.e. 0.001 GW cycles. 

This remarkably small dephasing concerns a "tuned" 
EOB phasing (with optimized flexibility parameters 05 
and i)poie)- However, as it is clear on Fig. 2, even the 
standard , "non-tuned" EOB phasing corresponding to 
our current analytical knowledge 05 = 0, z^poie = ''^piic(^)' 
agrees quite well with the NR phasing over the entire sim- 
ulation time. To exhibit this important fact in quantita- 
tive detail we compare in Fig. 6 the (splined) NR phase 
0nr(^') (after suitable shifts in (j) and t) to the stan- 
dard, "non-tuned" EOB phase (j)EOB{t) {0,5 — 0, Wpoic — 
''^poic(^))- ^^ ^^^'^ visual agreement (top panel) is too good 
to allow one to distinguish the two curves, we show (bot- 
tom panel) the phase difference (j)EOB{t) ~ 4>NR{t'). As 
expected, the dephasing is less good than in the above 
"tuned" case, but it remains impressively good: ±0.05 
radians, i.e., ±0.008 GW cycles, over the full time inter- 
val 900Af < tEOB < 3460M. 

Finally, we claim that, not only the phase, but also 
the amplitude of the new, resummed EOB waveform 



Eq. (4) exhibits a remarkable agreement with the NR 
data of [20]. Again, as Ref. [20] gave their results in 
the form of differences T4-NR, we plot in Fig. 7 the 
analog of the right panel of Fig. 19 there. We choose 
again the "optimum" values 05 — 40, Wpoic = 0.5074 
used in Fig. 5 and plot the NR -^ EOB analogs of the 
curves plotted by them in Fig. 19. Namely, we plot, 
at once, the W2-, W3- and W4-matched amplitude differ- 
ences [A'^"^/A]T4gQB = (A;^4 -^eob)Meob, where, 
as above, the T4 time is shifted so that uJT4:{t') and 
WEOB(i) agree when ujEOBitm) = ^m- In addition, 
we plot, as empty circles, some points taken (by ap- 
proximate measurements of ours) from the correspond- 
ing curve [A"*^/A]t4nr plotted on the right panel of 
Fig. 19 of [20]. The remarkable visual agreement be- 
tween these empty circles and our (A"* A/A)rp4gQB curve 
shows that: (i) the new, resummed 3 PN amplitude in- 
troduced in Ref. [30] and defined in Eqs. (4) (11) above 
agrees remarkably well with the NR one on the full 
time interval, 900M < iEOB < 3460M, (ii) by con- 
trast the Taylor T4 3.5/2.5 PN amplitude shows a signif- 
icant disagreement (~ —8%) in the same interval. Note 
that, though Ref. [20] emphasizes that the non-resummed 
3 PN-accurate waveform of [41] "improves agreement sig- 
nificantly" compared to the 2.5 PN one (used above), this 
improvement only concerns the early part of the inspiral. 
Indeed, Fig. 21 of [20] shows that the amplitude of Tay- 
lor T4 3.5/3.0 tends again to diverge together with Tay- 
lor T4 3.5/2.5 at the end of the inspiral: i.e., we think, 
precisely around the "dip" (near ^4) exhibited in Fig. 7 
above. 
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V. CONCLUSIONS 

We have investigated the agreement (in phase and 
in ampUtude) between the predictions of the EfFective- 
One-Body (EOB) formahsm and some accurate numer- 
ical data. We used as numerical data both (as a warm 
up) some old results on the energy flux from circular or- 
bits of a test mass around a non spinning black hole [39] , 
and some very recent results of the Caltech-Cornell group 
about the i = m = 2 gravitational wave emitted by 15 
orbits of an inspiralling system of two equal-mass non- 
spinning black holes [20] . 

In our warm up, test-mass example we showed how a 
slight tuning of the flexibility parameter [25] Wpoio (away 
from the naively expected value t'p*Jj;"'*''''*(i^ — 0) = 

1/V3 == 0.57735) to the value v^^fj:^ = 0) ~ 0.540 
allowed one to fit remarkably well the flux function 
F{v; v ~ 0) during the full inspiral, < v < wlso = 
1/V6. 

In the comparable mass case {ly — mim2/ {mi+m2)^ ^ 
1/4) we followed [30] in introducing a new, resummed 
3 PN-accuratC"'^^ EOB-type £ = m = 2 waveform. We 
then showed how to compute, for any values of the EOB 
flexibility parameters as (parametrizing 4 PN and higher 
conservative orbital interactions) and Vpoic (parametriz- 
ing i/-dependent 4 PN and higher effects in the re- 
summed radiation reaction) the EOB predictions for the 
£ = m — 2 gravitational curvature wave ^|^ ex d^h^^^ ex 
^EOB(t)e'''^^«^(*). 

We then compared the EOB predictions for the grav- 
itational wave (GW) phase, (/'EOB(i), and amplitude, 
^EOB(i), to the numerical relativity results of [20], say 
0NR(i), ^NR(i), using often as intermediary (as Ref. [20]) 
the so-called Taylor T4 3.5/2.5 post-Newtonian predic- 
tions (f>T4^(t), ATiit). Our main conclusions are: 

(i) In the GW frequency domain Muj < 0.08 where 
the Taylor T4 3.5/2.5 phase matches well with the NR 
phase, the EOB phase matches at least as well with the 
NR phase. A good EOB/NR match is obtained both 
for the standard "non-tuned" EOB flexibility parameters 
as = 0, Wpoio = ^piie('^) corresponding to our current an- 
alytical knowledge [1, 4] and for "tuned" EOB flexibility 
parameters. 

(ii) For higher GW frequencies, 0.08 < Moj < 0.14, 
while Taylor T4 3.5/2.5 starts to signiflcantly diverge 
from the NR phase, we showed that the standard "non- 
tuned" EOB phasing continues to stay in phase with NR 
within ±8 x 10^"^ GW cycles (see Fig. 6). Moreover, 
one can calibrate a^ and Wpoio so that the EOB phase 
matches with the NR phasing to the truly remarkable 
level of ±10-3 GW cycles over 30 GW cycles! 



^® Actually, our waveform has a greater accuracy than 3 PN in 
that it incorporates the test-mass limit of the 4 PN and 5 PN 
amplitude corrections. We shall occasionally refer to this PN 
accuracy as being S+'^-PN. 



(iii) We proposed several ways of "best fitting" the 
(flSj ^'poio)-dependent EOB predictions to accurate NR 
data: (a) by using the intrinsic representation of the 
phase evolution given by the reduced phase-acceleration 
function at^(w), Eq. (29); (b) by using selected ratios 
A'^"*(/)T4eob/A'^"*(/)T4nr and constraining them to be 
close to unity; and (c) by using an Lao norm of the dif- 
ference between (wm-matched) (f>EOB (t) and <j>'^'^ U'^ ) . 

Our results are given in several Figures. Notably, Fig. 3 
gives, for each given value of as, what is the optimum 
value of Wpoie which best fits (in the sense of the ratio 
/9^7<^, Eq. (35)) the NR data. Then, Fig. 4 plots the Loo 
distance (on a large time-interval roughly correspond- 
ing to the full simulation of [20]) betweem (/)EOB(i) and 
'/'nr (C4) as a function of as (for Wpoio = f polo (as) given 
by Fig. 3). We find that the absolute value of the max- 
imum dephasing between EOB and NR can be as small 
as 0.006 radians (or 0.001 GW cycles) if 30 < as < 52. 
However, it is difficult to be precise about the "preferred" 
valued of as. We recall in this respect that, recently, 
Ref. [29] has tried to constrain the value of as (keep- 
ing, however, Upoie fixed to v^lf^{i'), and without using 
our improved EOB waveform) by maximizing the over- 
lap between EOB and NR plunge waveforms. They found 
that the overlap was good (and fiat) over a rather large 
interval of values of as (that they denote as A), roughly 
centered around as — 60. We note, however, that this be- 
havior might be due (at least in part) to the phenomenon 
pointed out in [25]. In the latter reference (where as was 
denoted as 65), it was found that the use of EOB tem- 
plates based on as = 50 (rather than as — 0) allowed one 
to have large overlaps (large "effectualnesses" ) with all 
other EOB templates. At this stage, we therefore do not 
have yet any precise knowledge of what might be the pre- 
ferred "effective" value of as. Our work, however, shows 
that there is a quite strict correlation between the best-fit 
choices of as and fpoic- When, in the future, 05 becomes 
precisely known, it will be interesting to see what is the 
corresponding value of Upoie('^ = 1/4) and to compare it 
to the best-fit value Wpoio(i^ = 0) ~ 0.540 obtained in our 
warm-up Sec. II. 

For instance, the couple as = 40, Wpoic = 0.5074 yields 
a remarkable good fit to the NR data reported in [20]. 
We show the comparison of the various phasings (NR, 
EOB, T4) in Fig. 5. This Figure clearly exhibits how our 
best-fit EOB phase does a much better job than any non- 
resummed PN approximant at following the NR phase. 
We finally get dephasings smaller than ±0.006 radians 
(i.e. < 10~3 GW cycles!) over about 30 GW cycles! 

Finally, we exhibited in Fig. 7 how the amplitude of 
our new, resummed 3"'"^-PN-accurate EOB waveform, 
Eq. (4), exhibits a remarkable agreement with the corre- 
sponding amplitude of the NR data of [20]. The agree- 
ment is clearly better than any, non resummed PN am- 
plitude, including the recent 3 PN-accurate one of Kid- 
der [41]. 

We think that the present work, taken in conjunc- 
tion with other recent works on the EOB-NR compar- 
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ison [29] [28, 30], confirms the remarkable ability of 
the EOB formalism (especially in its recently improved 
avatars) to agree with NR results. Note in particular 
that the level of phase agreement reached here is bet- 
ter by a factor 30 (±0.001 GW cycles versus ±0.03 GW 
cycles for v = 1/4) than what was recently achieved, 
for merger signals, in Ref. [29] using less accurate ver- 
sions of EOB waveforms than the one used here. We 
suggest that the ground-based interferometric GW de- 
tectors should include in their template banks the new, 
extended and improved EOB waveforms which are be- 
ing developed and notably the resummed one introduced 
in [30] and generalized here. We also suggest that NR 
data be made available in some repository, soon after 



the first published results, to expert theorists willing to 
extract the physical information they contain. 
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